A biallelic multiple nucleotide length polymorphism explains functional causality at 5p15.33 prostate cancer risk locus.

To date, single-nucleotide polymorphisms (SNPs) have been the most intensively investigated class of polymorphisms in genome wide associations studies (GWAS), however, other classes such as insertion-deletion or multiple nucleotide length polymorphism (MNLPs) may also confer disease risk. Multiple reports have shown that the 5p15.33 prostate cancer risk region is a particularly strong expression quantitative trait locus (eQTL) for Iroquois Homeobox 4 (IRX4) transcripts. Here, we demonstrate using epigenome and genome editing that a biallelic (21 and 47 base pairs (bp)) MNLP is the causal variant regulating IRX4 transcript levels. In LNCaP prostate cancer cells (homozygous for the 21 bp short allele), a single copy knock-in of the 47 bp long allele potently alters the chromatin state, enabling de novo functional binding of the androgen receptor (AR) associated with increased chromatin accessibility, Histone 3 lysine 27 acetylation (H3K27ac), and ~3-fold upregulation of IRX4 expression. We further show that an MNLP is amongst the strongest candidate susceptibility variants at two additional prostate cancer risk loci. We estimated that at least 5% of prostate cancer risk loci could be explained by functional non-SNP causal variants, which may have broader implications for other cancers GWAS. More generally, our results underscore the importance of investigating other classes of inherited variation as causal mediators of human traits.

Rewiring of master transcription factor cistromes during high-grade serous ovarian cancer development.

The transcription factors MECOM, PAX8, SOX17 and WT1 are candidate master regulators of high-grade serous 'ovarian' cancer (HGSC), yet their cooperative role in the hypothesized tissue of origin, the fallopian tube secretory epithelium (FTSEC) is unknown. We generated 26 epigenome (CUT&TAG, CUT&RUN, ATAC-seq and HiC) data sets and 24 profiles of RNA-seq transcription factor knock-down followed by RNA sequencing in FTSEC and HGSC models to define binding sites and gene sets regulated by these factors in cis and trans. This revealed that MECOM, PAX8, SOX17 and WT1 are lineage-enriched, super-enhancer associated master regulators whose cooperative DNA-binding patterns and target genes are re-wired during tumor development. All four TFs were indispensable for HGSC clonogenicity and survival but only depletion of PAX8 and WT1 impaired FTSEC cell survival. These four TFs were pharmacologically inhibited by transcriptional inhibitors only in HGSCs but not in FTSECs. Collectively, our data highlights that tumor-specific epigenetic remodeling is tightly related to MECOM, PAX8, SOX17 and WT1 activity and these transcription factors are targetable in a tumor-specific manner through transcriptional inhibitors.

Single-cell transcriptomic analysis of endometriosis.

Endometriosis is a common condition in women that causes chronic pain and infertility and is associated with an elevated risk of ovarian cancer. We profiled transcriptomes of >370,000 individual cells from endometriomas (n = 8), endometriosis (n = 28), eutopic endometrium (n = 10), unaffected ovary (n = 4) and endometriosis-free peritoneum (n = 4), generating a cellular atlas of endometrial-type epithelial cells, stromal cells and microenvironmental cell populations across tissue sites. Cellular and molecular signatures of endometrial-type epithelium and stroma differed across tissue types, suggesting a role for cellular restructuring and transcriptional reprogramming in the disease. Epithelium, stroma and proximal mesothelial cells of endometriomas showed dysregulation of pro-inflammatory pathways and upregulation of complement proteins. Somatic ARID1A mutation in epithelial cells was associated with upregulation of pro-angiogenic and pro-lymphangiogenic factors and remodeling of the endothelial cell compartment, with enrichment of lymphatic endothelial cells. Finally, signatures of ciliated epithelial cells were enriched in ovarian cancers, reinforcing epidemiologic associations between these two diseases.

Predicting master transcription factors from pan-cancer expression data.

Critical developmental "master transcription factors" (MTFs) can be subverted during tumorigenesis to control oncogenic transcriptional programs. Current approaches to identifying MTFs rely on ChIP-seq data, which is unavailable for many cancers. We developed the CaCTS (Cancer Core Transcription factor Specificity) algorithm to prioritize candidate MTFs using pan-cancer RNA sequencing data. CaCTS identified candidate MTFs across 34 tumor types and 140 subtypes including predictions for cancer types/subtypes for which MTFs are unknown, including e.g. PAX8, SOX17, and MECOM as candidates in ovarian cancer (OvCa). In OvCa cells, consistent with known MTF properties, these factors are required for viability, lie proximal to superenhancers, co-occupy regulatory elements globally, co-bind loci encoding OvCa biomarkers, and are sensitive to pharmacologic inhibition of transcription. Our predictions of MTFs, especially for tumor types with limited understanding of transcriptional drivers, pave the way to therapeutic targeting of MTFs in a broad spectrum of cancers.

Single-cell transcriptomics identifies gene expression networks driving differentiation and tumorigenesis in the human fallopian tube.

The human fallopian tube harbors the cell of origin for the majority of high-grade serous "ovarian" cancers (HGSCs), but its cellular composition, particularly the epithelial component, is poorly characterized. We perform single-cell transcriptomic profiling of around 53,000 individual cells from 12 primary fallopian specimens to map their major cell types. We identify 10 epithelial subpopulations with diverse transcriptional programs. Based on transcriptional signatures, we reconstruct a trajectory whereby secretory cells differentiate into ciliated cells via a RUNX3high intermediate. Computational deconvolution of advanced HGSCs identifies the "early secretory" population as a likely precursor state for the majority of HGSCs. Its signature comprises both epithelial and mesenchymal features and is enriched in mesenchymal-type HGSCs (p = 6.7 × 10-27), a group known to have particularly poor prognoses. This cellular and molecular compendium of the human fallopian tube in cancer-free women is expected to advance our understanding of the earliest stages of fallopian epithelial neoplasia.

A Study of High-Grade Serous Ovarian Cancer Origins Implicates the SOX18 Transcription Factor in Tumor Development.

Fallopian tube secretory epithelial cells (FTSECs) are likely the main precursor cell type of high-grade serous ovarian cancers (HGSOCs), but these tumors may also arise from ovarian surface epithelial cells (OSECs). We profiled global landscapes of gene expression and active chromatin to characterize molecular similarities between OSECs (n = 114), FTSECs (n = 74), and HGSOCs (n = 394). A one-class machine learning algorithm predicts that most HGSOCs derive from FTSECs, with particularly high FTSEC scores in mesenchymal-type HGSOCs (padj < 8 × 10-4). However, a subset of HGSOCs likely derive from OSECs, particularly HGSOCs of the proliferative type (padj < 2 × 10-4), suggesting a dualistic model for HGSOC origins. Super-enhancer (SE) landscapes were also more similar between FTSECs and HGSOCs than between OSECs and HGSOCs (p < 2.2 × 10-16). The SOX18 transcription factor (TF) coincided with a HGSOC-specific SE, and ectopic overexpression of SOX18 in FTSECs caused epithelial-to-mesenchymal transition, indicating that SOX18 plays a role in establishing the mesenchymal signature of fallopian-derived HGSOCs.

Expression of the BAD pathway is a marker of triple-negative status and poor outcome.

Triple-negative breast cancer (TNBC) has few therapeutic targets, making nonspecific chemotherapy the main treatment. Therapies enhancing cancer cell sensitivity to cytotoxic agents could significantly improve patient outcomes. A BCL2-associated agonist of cell death (BAD) pathway gene expression signature (BPGES) was derived using principal component analysis (PCA) and evaluated for associations with the TNBC phenotype and clinical outcomes. Immunohistochemistry was used to determine the relative expression levels of phospho-BAD isoforms in tumour samples. Cell survival assays evaluated the effects of BAD pathway inhibition on chemo-sensitivity. BPGES score was associated with TNBC status and overall survival (OS) in breast cancer samples of the Moffitt Total Cancer Care dataset and The Cancer Genome Atlas (TCGA). TNBC tumours were enriched for the expression of phospho-BAD isoforms. Further, the BPGES was associated with TNBC status in breast cancer cell lines of the Cancer Cell Line Encyclopedia (CCLE). Targeted inhibition of kinases known to phosphorylate BAD protein resulted in increased sensitivity to platinum agents in TNBC cell lines compared to non-TNBC cell lines. The BAD pathway is associated with triple-negative status and OS. TNBC tumours were enriched for the expression of phosphorylated BAD protein compared to non-TNBC tumours. These findings suggest that the BAD pathway it is an important determinant of TNBC clinical outcomes.

Sensitivity of ovarian cancer cells to acetaminophen reveals biological pathways that affect patient survival.

Experimental and epidemiological data support the potential activity of acetaminophen against ovarian cancer (OVCA). In this study, we sought to confirm the activity of acetaminophen in OVCA cell lines and to investigate the molecular basis of response. A total of 16 OVCA cell lines underwent pretreatment (baseline) genome-wide expression measurements and were then treated with and analyzed for acetaminophen sensitivity. Pearson's correlation analysis was performed to identify genes that were associated with OVCA acetaminophen response. The identified genes were subjected to pathway analysis, and the expression of each represented pathway was summarized using principal component analysis. OVCA acetaminophen response pathways were analyzed in 4 external clinico-genomic datasets from 820 women for associations with overall survival from OVCA. Acetaminophen exhibited antiproliferative activity against all tested OVCA cell lines, with half maximal inhibitory concentration values ranging from 63.2 to 403 µM. Pearson's correlation followed by biological pathway analysis identified 13 pathways to be associated with acetaminophen sensitivity (P<0.01). Associations were observed between patient survival from OVCA and expression of the following pathways: Development/angiotensin signaling via ß-arrestin (P=0.04), protein folding and maturation/angiotensin system maturation (P=0.02), signal transduction/c-Jun N-terminal kinase (JNK) pathway (P=0.03) and androstenedione and testosterone biosynthesis and metabolism (P=0.02). We confirmed that acetaminophen was active against OVCA cells in vitro. Furthermore, we identified 4 molecular signaling pathways associated with acetaminophen response that may also affect overall survival in women with OVCA, including the JNK pathway, which has been previously implicated in the mechanism of action of acetaminophen and is predictive of decreased survival in women with OVCA.

Micro-RNAs associated with the evolution of ovarian cancer cisplatin resistance.

OBJECTIVES: Ovarian cancer (OVCA) is the leading cause of mortality among women with gynecologic malignancy, in part due to the development of chemoresistance. We sought to identify micro-RNAs (miRNAs) associated with in vitro development of OVCA chemoresistance that may also represent potential targets for therapy. METHODS: In this study, four OVCA cell lines (A2780CP, A2780S, IGROV1, and OVCAR5) were serially treated with cisplatin in parallel with measurements of miRNA expression changes. RESULTS: Nine miRNAs were found to be associated with increasing cisplatin resistance (IC50) (p<0.01); however, only 5 of these miRNAs have publically available information. Pathway analysis identified 15 molecular signaling pathways that were represented by genes predicted to be targets of the 5 miRNAs (false discovery rate<0.05), 11 of which are associated with the epithelial-mesenchymal transition (EMT). Further analysis identified 2 of those pathways as being associated with overall survival in 218 patients with OVCA. CONCLUSIONS: Collectively, this panel of miRNAs associated with in vitro evolution of OVCA cisplatin resistance and the pathways identified to be associated with EMT and overall patient survival provide a framework for further investigations into EMT as a therapeutic target in patients with OVCA.

BAD-mediated apoptotic pathway is associated with human cancer development.

The malignant transformation of normal cells is caused in part by aberrant gene expression disrupting the regulation of cell proliferation, apoptosis, senescence and DNA repair. Evidence suggests that the Bcl-2 antagonist of cell death (BAD)-mediated apoptotic pathway influences cancer chemoresistance. In the present study, we explored the role of the BAD-mediated apoptotic pathway in the development and progression of cancer. Using principal component analysis to derive a numeric score representing pathway expression, we evaluated clinico-genomic datasets (n=427) from corresponding normal, pre-invasive and invasive cancers of different types, such as ovarian, endometrial, breast and colon cancers in order to determine the associations between the BAD-mediated apoptotic pathway and cancer development. Immunofluorescence was used to compare the expression levels of phosphorylated BAD [pBAD (serine-112, -136 and -155)] in immortalized normal and invasive ovarian, colon and breast cancer cells. The expression of the BAD-mediated apoptotic pathway phosphatase, PP2C, was evaluated by RT-qPCR in the normal and ovarian cancer tissue samples. The growth-promoting effects of pBAD protein levels in the immortalized normal and cancer cells were assessed using siRNA depletion experiments with MTS assays. The expression of the BAD-mediated apoptotic pathway was associated with the development and/or progression of ovarian (n=106, p<0.001), breast (n=185, p<0.0008; n=61, p=0.04), colon (n=22, p<0.001) and endometrial (n=33, p<0.001) cancers, as well as with ovarian endometriosis (n=20, p<0.001). Higher pBAD protein levels were observed in the cancer cells compared to the immortalized normal cells, whereas PP2C gene expression was lower in the cancer compared to the ovarian tumor tissue samples (n=76, p<0.001). The increased pBAD protein levels after the depletion of PP2C conferred a growth advantage to the immortalized normal and cancer cells. The BAD-mediated apoptotic pathway is thus associated with the development of human cancers likely influenced by the protein levels of pBAD.

The Chinese herb polyphyllin D sensitizes ovarian cancer cells to cisplatin-induced growth arrest.

PURPOSE: We evaluated the effects of polyphyllin D (PD), a natural compound with anti-neoplastic activity and a major component of the Chinese herb Paris polyphylla, on ovarian cancer (OVCA) cell line proliferation and platinum sensitivity. METHODS: A panel of 20 OVCA cell lines was subjected to PD treatment, MTS proliferation assays, and determination of IC50. Pre-treatment, baseline genome-wide Affymetrix expression analysis was performed on each cell line, and Pearson's correlation was performed to identify genes associated with OVCA PD sensitivity. Twelve cell lines were treated with PD with and without cisplatin, and the effects of PD on cisplatin IC50 were quantified. Genes associated with OVCA PD sensitivity were evaluated for associations with survival in a publically available clinico-genomic dataset of 218 patients with OVCA. RESULTS: Our results showed that PD exhibited anti-proliferative effects against all OVCA cell lines tested, with IC50 values ranging from 0.2 to 1.4 µm. Furthermore, in all cell lines, PD treatment significantly decreased cisplatin IC50 (mean IC50 reduction of 2.1 µm; P < 0.02). Pearson's correlation test identified 25 probe sets, representing 18 unique genes to be associated with PD sensitivity (FDR = 0). We found that one of these genes was associated with overall survival in women with OVCA: CLDN4 (P = 0.014). CONCLUSION: Our findings highlight the value of PD as a natural product with anti-cancer properties, which may also enhance the activity of existing therapeutic agents.

Gene expression data reveal common pathways that characterize the unifocal nature of ovarian cancer.

OBJECTIVE: The objective of the study was to evaluate the biological validity of ovarian cancer (OVCA) screening and early detection efforts and to characterize signaling pathways associated with human cancer metastasis and patient survival. STUDY DESIGN: Using genome-wide expression profiling and deoxyribonucleic acid sequencing, we compared pelvic and matched extrapelvic implants from 30 patients with advanced-stage OVCA for expression of molecular signaling pathways and p53 gene mutations. Differentially expressed pathways were further evaluated in a series of primary or early-stage vs metastatic or recurrent cancer samples from 389 ovarian, prostate, and oral cancer patients. Metastasis pathways were also evaluated for associations with survival in 9 independent clinicogenomic datasets from 1691 ovarian, breast, colon, brain, and lung cancer and leukemia patients. The inhibitory effects of 1 pathway (transforming growth factor [TGF]-WNT) on in vitro OVCA cell migration were studied. RESULTS: Pelvic and extrapelvic OVCA implants demonstrated similar patterns of signaling pathway expression and identical p53 mutations. However, we identified 3 molecular pathways/cellular processes that were differentially expressed between pelvic and extrapelvic OVCA samples and between primary/early-stage and metastatic/advanced or recurrent ovarian, oral, and prostate cancers. Furthermore, their expression was associated with overall survival from ovarian cancer (P = .006), colon cancer (1 pathway at P = .005), and leukemia (P = .05). Artesunate-induced TGF-WNT pathway inhibition impaired OVCA cell migration. CONCLUSION: Advanced-stage OVCA has a unifocal origin in the pelvis. Molecular pathways associated with extrapelvic OVCA spread are also associated with metastasis from other human cancers and with overall patient survival. Such pathways represent appealing therapeutic targets for patients with metastatic disease.

Perifosine, an AKT inhibitor, modulates ovarian cancer cell line sensitivity to cisplatin-induced growth arrest.

OBJECTIVES: AKT, a key regulator of diverse tumor signaling, is associated with progression of many cancers. Here, we investigated 1) the influence of AKT on survival from ovarian cancer (OVCA), 2) the activity of the AKT inhibitor perifosine ± cisplatin, and 3) the molecular determinants of perifosine-response. Phospho-AKT expression values and Affymetrix U133a expression data were downloaded from The Cancer Genome Atlas. METHODS: Pearson correlation was used to determine associations between overall survival from OVCA and therapy response. Genes and represented signaling pathways associated with perifosine-response were explored in OVCA cells (n=10) and the NCI60 cancer cell panel. Pathway expressions, modeled by PCA, were evaluated for influences on survival using publically available clinico-genomic datasets. RESULTS: Phospho-AKT (serine473) expression correlated with survival from OVCA (P<0.05) and platinum-response (P=0.004). In vitro, perifosine showed anti-proliferative effects against OVCA cells and potentiated cisplatin-induced growth arrest. Perifosine-response was associated with the expression (FDR<0.05) of 7 signaling pathways in OVCA cells and 64 signaling pathways in the NCI60 cell panel. Three pathways were found in common: 1) Cytoskeleton remodeling/cytoskeleton remodeling (cyto), 2) cell adhesion/chemokines and adhesion (chemokines), and 3) cytoskeleton remodeling/TGF-WNT (TGF-WNT). The TGF-WNT was associated with survival from OVCA (P=0.0055). CONCLUSIONS: AKT signaling is an important determinant of OVCA response to chemotherapy and overall patient survival. Our data provide insight into the molecular basis to perifosine activity and identifies pathways associated with perifosine sensitivity and patient clinical outcome.